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ABSTRACT 

This work examines the spin-axis attitude stability of dissipative dual-spin 
spacecrafts in force-free environment. A mathematical model of a representa- 
tive dual-spin system with energy dissipation on both bodies is proposed and the 
dynamical equations of motion are derived. The linear variational equations 
describing the behavior of the perturbed system in a neighborhood of an assumed 
equilibrium solution are shown to have periodic coefficients. A suitable similar- 
ity transformation is introduced which reduces (in the sense of Liapunov) the 
time-varying system to an autonomous one. This allows all the powerful tech- 
niques of time-invariant linear system theory to be applied to the dual-spin sys- 
tem. The theory is applied on a recently developed mathematical model of the 
SAS-A spacecraft and it is demonstrated that this spacecraft may be analyzed 
more precisely by utilizing this transformation than what has so far been 
achieved. 
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LIST OF SYMBOLS 

0123: reference frame aligned with the principal axes of inertia 
A(t): system matrix of an arbitrary time-varying linear system 
afa^: moment arms for dampers on main body (rotor) 

B: transformation matrix 

c(c'): viscous friction constant for dampers on main body (rotor) 
c (c 1 ): c(c* ) per unit mass 

d(d'): axis 3 location of centers of mass of unperturbed dampers on main 
body (rotor) 

h: angular momentum 

I: matrix representation of the inertia tensor of the system, 

1 = (fjic); j = i. 2 . 3 ; k= !» 2. 3 

Ii , I 2 , 1 3 : principal moments of inertia of the system 4*- 

J 3 (J' ): spin axis moment of inertia of the main body (rotor) 

i 

k(k'): spring constant for dampers on the main body (rotor) . 

k(k'): k(k') per unit mass 

l 

■i (l ' ): axis 3 location of center of mass of the main body (rotor) 

Lj, L 2 (Lj, L 2 ): system parameters defined by L, = ma/l,, 

L 2 = ma/lj, Lj = m'a'/li, L 2 = m t a. , /l 2 j 

M x : total mass of the system, M T = M + M' + 4m + 4m’ 

1 
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M(M'); mass of the main body (rotor) 

m, , m 2 (mj , m'): masses of dampers on the main body (rotor) 
m,, m. (m* , m'): balancing masses on the main body (rotor) 



m(m' ): m - = m 2 = m 3 = m 4 (m* = m| = m' 2 = m' 3 = m' 4 ) 

n: nominal spin rate of the main body about axis 3 

nj : n t = n + q 

Uj : n 2 = n + 2q 

P(t): transformation matrix 

Q(s): system matrix in the complex frequency domain 
s: complex frequency, Laplace transform variable 

X: varying part of transverse axes moments of inertia during perturbed motion 

y: axis 3 motion of the main body and the rotor during perturbed motion (with 
respect to the center of mass of the system) 

y r , y 2 : transformed variables related to the displacements of the rotor dampers 

z i * z 2 ( z i* Z 2 displacements, with respect to null positions, of dampers on the 
main body (rotor) 

k j , k 2 : frequencies 

, a> 2 , cc 3 : system angular rates 

: skew symmetric matrix involving angular rates 

p(p'): mass ratio, p = m/M T ( p ' = m’/M T ) 

6 : angular orientation of the rotor, about axis 3, with respect to the main body 
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INTRODUCTION 


It is well known that an ideal rigid body rotating in force -free environment 
is stable in spin about the axis of either the greatest or the least principal moment 
of inertia. Where spacecraft applications permit, spin stabilization is an extremely 
attractive and widely used technique for passive attitude stabilization. Following 
the launch and subsequent unstable behavior demonstrated by Explorer I in 1958 
came the realization that internal energy dissipation may cause the spinning 
motion about the axis of the least principal moment of inertia to become unstable. 
To get around this problem spinning spacecrafts are often designed to spin about 
their axes of the greatest principal moments of inertia. Where this cannot be 
done for some reason, the so called "dual-spin" configuration may be employed. 

The basic idea is to attach a relatively large spinning flywheel rigidly to the 
main spacecraft which may spin, if at all, at a much slower rate. Since the total 
angular momentum remains inertially fixed in the assumed force-free environ- 
ment, it is expected that the system spin axis would remain attitude stabilized 
if the main spacecraft momentum is small. 

For the lossless system, it is a simple matter to extend the classical anal- 
ysis for the Liapunov stability of a single rigid spinning body to the dual-spin 
case. The main question of interest is of course the stability of the system 
with the energy dissipation taken into account. An extensive literature exists 
reporting the results of various design studies of specific dual-spin systems 
and damper configurations. Not much analytical work has however been done 

1 


to explain satisfactorily the role of energy dissipation in causing instability of 
steady spinning motion. Likins 1 was one of the first to derive a quantitative 
trade-off relationship between the stabilizing and destabilizing effects of energy 
dissipations in the two component parts of a dual-spin system. He applied a 
heuristic energy sink approximation to relate momentum variation with energy 
dissipation: a very valuable result from a spacecraft designer’s point of view. 
However, he emphatically disclaims rigor and Mingori 2 states that detailed 
numerical analysis shows that the results obtained from energy sink approxi- 
mations can in certain cases be obscure or even misleading. 


Likins and Mingori 3 discuss some of the considerations involved in a proper 
choice of variables. Since the inertial attitude stability of the spin axis is of 
interest, the natural choice is a set of attitude coordinates (e.g., Euler angles) 
in an inertial frame attached to the center of mass of the system. For the gen- 
eral case cf a mechanical system with several connected parts capable of storing 
angular momenta in various orientations the equations must indeed be formulated 
in terms of attitude coordinates. Pringle 4 has done this to develop a general 
stability theory of mechanical systems but it is not quite clear if his theory can 
be applied to a dual-spin system with energy dissipation on both bodies. 

The equations of motion are more conveniently written in terms of the angular 
rates resolved along a body fixed reference frame. It can be easily seen that for the 
specific case of the dual-spin configuration, some form of Euler's equations in 
terms of angular rates suitably augmented by energy dissipation aquations is an 


2 

•rtinTHh HPttir fraa n r w a r » » • - -- 




1 

I 

i 


j i 



adequate mathematical model. Ir this case, the equilibrium motion is a spin 
about one of the body axes which is the common axis of spin of all the connected 
parts of the system. In stable spin equilibrium the total system angular momentum 
is fixed in the body and inertial attitude stability of the spin axis may be inferred 
by an appeal to the fact that this momentum vector is inertially fixed in the 
assumed force-free environment. 

Accurate modeling is extremely difficult to achieve for the slightly flexible 
dissipative spacecraft structures and some of the various damping mechanisms 
that have been proposed or used in dual-spin systems. Precise stability results 
can however be obtained only for mathematically well defined problems. We 
therefore look for a model which is of some practical interest in itself and is 
also general enough to serve as an approximate model for a large class of dual- 
spin systems. Mingori 2 has considered such a model and we take a somewhat 
more general version of his system as our basic dissipative dual-spin config- 
uration. In this model, the energy dissipation in any component part, whether 
inherent or introduced through a damping mechanism, is modeled as produced 
by one or more ball-in-tube dampers. This damper is a practical realization 
of an ideal linear damped harmonic oscillator and consists of a spring-restrained 
point mass in a short straight viscous fluid filled tube. The generalization from 
Mingori's model lies in the appreciation of the fact that the high speed rotor is 
always made symmetric and the rotor energy dissipation must therefore be 
symmetric in the two transverse axes. 



The equations o f motion of the basic dissipative dual-spin system are highly 
nonlinear. Rigorous approaches, such as an attempt to construct a Liapunov func- 
tion, do not appear to be very promising. The natural attack on the problem is to 
apply the Principle of Stability in the First Approximation 5 which succeeds easily 
if the linearized variational system of equations of perturbed motion is time- 
invariant and does not have purely Imaginary characteristic roots. It turns out 
however that the linearized equations for this system have periodic coefficients. 
As Hahn 5 points out, the known stability results for such systems are not exactly 
stability criteria because in order to apply them several things must already be 
known about the solutions. 

An exact method of analysis may be developed by use of Floquet's theorem 
and the Liapunov Reducibility theorem . 5 Each linear differential equation with 
periodic coefficients can be reduced (in the sense of Liapunov) to a kinematically 
similar autonomous differential equation having the same stability behavior by 
use of a nonsingular matrix with periodic coefficients. Stability and performance 
analysis of the reduced system is a fairly routine exercise in linear system 
theory. The main difficulty in applying this approach is of course finding the 
appropriate similarity transformation. For the proposed dissipative dual-spin 
model physical Insight suggests a very simple quasi -holonomic transformation 
which serves the purpose. 
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MATHEMATICAL MODEL OF ^HE BASIC DISSIPATIVE 
DUAL-SPIN CONFIGURATION 

Consider the model in Fig. 1 which shows an asymmetrical rigid body of 
mass M, which will be called the main body, carrying a uniaxial rotor of mass 
M* sp inning about a principal axis. Assume that the energy dissipation on the 
rotor can be modelled by two ball-in-tube dampers of mass mj and mj , spring 
constant c* and viscous friction coefficient k\ each with a single degree of 
freedom in a direction parallel to the wheel axis. I jet there be two other mass 
points m' 3 and mj for static and dynamic balance of the system around the spin 
axis. To counteract the destabilizing effect of rotor dissipation, let there be 
two similar dampers with masses m, and m 2 , spring constant c and viscous 
friction coefficient k, as well as two balancing masses m 3 and m 4 on the main 
body. The configuration is more precisely described by the mathematical 
formulation below. 

Let 0 be the center of mass of the system, which is fixed in inertial space, 
and consider a rectangular frame 0123 along the principal axes of inertia of the 
system at rest. Let 3 be the spin axis. Let 6 be the angular orientation of the 
rotor relative to the main body and assume that an ideal servo maintains 8 con- 
stant at a value q. Hie table below specifies the centers of masses c ' the vari- 
ous bodies both when the system is at rest and when it is undergoing perturbed 
motion. 
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Table 1. Configuration of the Basic Dual-Spin System 


Mass Rest Position Perturbed Position 


M 

0 

0 

i 

0 

0 

-£ + y 

M' 

0 

0 

V 

0 

0 

V + y 

m ! 

a 

0 

d 

a 

0 

d + y + Zj 

m 2 

0 

a 

d 

0 

a 

d + y + z 2 

m 3 

-a 

0 

d 

-a 

0 

d + y 

m 4 

0 

-a 

d 

0 

-a 

d + y 

m l 

a' 

0 

d' 

a 1 cos 6 

a' sin 0 

d + y + zj 

m 2 

0 

a’ 

d' 

-a 1 sin 8 

a 1 cos# 

d' : y + z' 

m 3 

-a' 

0 

d' 

-a # cos 6 

-a # sin# 

d' + y 

m 4 

0 

-a' 

d' 

a' sin# 

-a' cos# 

d' + y 


For simplicity, we have assumed that m t = m 2 = m 3 = m 4 - m and mj = m' 2 
= m\ -= m' = m'. The damper mass displacements from their respective null 

w 4 

positions are z,, z 2 , z\ and z\. Equating moments about the origin 

y = -p(z j + z 2 ) - P* (z' + z' ) (1) 

where p = m/M T ,p'= m'/M T and = M + M’ + 4m + 4m’ is the total mass of 
the system. 
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Let I jt I 2 and I 3 be the principal moments of inertia of the system at rest. 
Write I = J, + J’ where J, and J* are the moments of inertia of the main bodv 

3 3 3 3 3 

and the rotor, each including the corresponding dampers and balancing masses 
at rest. During perturbed motion, the components of the inertia tensor are 


I 


u 


= I. 


+ X 


i 22 = i 2 + x 

^33 = h = h + j ; 


1 23 .= I 32 = -ma z 2 - tn a 1 (z^ s in 6 + z ' 2 cos 0 ) 
I 13 = I 31 = -ma Zj - m* a' (zj cos 9 - z ' 2 sin#) 


where 


X = Mj.y 2 + m{z 2 + z 2 + 2 (z, + z 2 ) (d + y)} 

+ m* {zj 2 + z' 2 + 2 (z * + z') (d' + y)} 


( 2 ) 


The angular momentum of the system with the rotor and the dampers 
frozen in is given by I« where i is the matrix representation oi the inertia 
tensor defined above and « is the angular velocity vector, a?= col (co lt a> 2 , w 3 ). 
The rotor contributes an additional momentum of Jj q along the third axis. 
The calculation of the inertif 1 1 already takes into account the effects of much 
of the linear momenta of the damper masses, the additional terms to be 
considered are due to the time variations of the damper displacements and 
wheel rotation. From basic mechanics, a unit mass point at position 


x = col (Xj , x 2 , x 3 ) on a body with angular velocity a> = col (w lt w, , oj 3 ) contributes 
to the total angular momentum 


Ah = 


| x 2 £ 3 

~ X 3 X 2 

+ "l ( X 2 

+ X 2 ) - W 2 XjX 2 

“ w 3 x l X 3 

X 3 X 1 

* X 1 X 3 

-"l X 2 X l 

+ u 2 (x 2 + X 2 ) 

- W 3 X 2 X 3 

X 1 X 2 

" X 2 X 1 

-V3 X 1 

- <o 2 x 3 x 2 + w 3 

(x 2 + X 2 ) 


( 3 ) 


We sum the above expression over the dampers neglecting the terms already 
taken into account and add to I&> to get the total angular momentum of the system 


as 


h = 



+ X CcJ J 

+ ma (z 2 - oj 3 Zj) 




+ m # a 1 

’ {cos # (zj - (a; 3 + q) z! 

[) + s in t? (z' 

+ ( w 3 + q) zj)} 

I 2 cj 2 

4- Xu 2 

- ma (Zj 4- ct> 3 z 2 ) 




9 9 

- m a 

{-sin# (z 2 - (oj 3 + q) 

z|) + cos 8 (z 

; + (w 3 + q)zj)} 

I 3 Cd 3 

+ j;q 

- ma (co | z j + a > 2 z 2 / 




-m a* 

{cos #(,,. 7 ; +r., 2 z') - 

s i n 6 (gj z' “ 

Si 

“'2*1 7 


( 4 ) 


Inertial differentiation of the angular momentum yields the rotational equa- 
tions of motion as 

*1 "l ” {(*2 " * 3 ) ^3 * J 3 q } + X ("1 “ w 2 " 3 ) + * (0 l 

{z 2 - (fl) $ + «j6> 2 ) Zj - (a. 2 - cj\) z 2 } 
n'a' jcos#[zj - (<i 3 +a,w 2 ) z' - {<o 2 - («u 3 + q) 2 } z'J 
in# [z'j + (<i> 3 + co l oj 2 ) z' 2 - {w 2 - (w 3 + q) 2 } z{J| = 0 


■f ma 
+ m' 


+ s in 1 


( 5 ) 


1 

l 
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jg ****** 


I 2^2 * {( X 3 " X l> u 3 + Ja a 1 1 X (" 2 + W 3 w l) + Xw 2 

-ma (zj + (a> 3 - « t o> 2 ) z 2 + (a> 2 - <y 2 ) 2 J 
-mV | cos 6 [z; + (o>3 - Cv x m 2 ) z ' 2 .[(0)3 + q) 2 - 0J 2 }zjJ 

- sin^ £z' - (o> 3 - WjW 2 ) z[ + {(^’3 + q) 2 -&>i} z'jj = 0 (6) 

I 3^3 * < X i " X 2> 

-ma {2cdj z x + 2co 2 z 2 + (o>j - &> 2 ;a 3 ) 2 i + (^2 + w j w i) *2} 

-mV |cos 0 [2 co^j + 2 o> 2 Zj + (i, - 0^0)3) zj + (a> 2 + WjWj) z 2 ] 

- sin 8 [2o)j z' - 2<a 2 z' + (oj^ - u> 2 a> 3 ) z' - (Z > 2 + a^ajj) zj ]| = 0 00 


These equations have to be augmented by four others for the four dampers. 
From the basic rule for inertial differentiation, 6 a mass point with position 
x = col(Xj, x 2 , X3) has a linear acceleration in inertial space given by 

x +2fix+(fl+ fl 2 ) x (8) 


where 


- 0 ), 


n = 


o>. 


( 9 ) 


The equation of motion for any damper is obtained by summing the third com- 
ponent of the vector expression (8) with the corresponding viscous and elastic 
forces and equating to zero. Thus 

^y' + z 1 -a(w 2 - 6 , 3 « 1 )-(^ 2 +a ) 2 ) (d + y + Zj)} H-cZj+kZjrO ( 10 ) 


ml 
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m { y + z 2 + a (cl) x + a; 2 6; 3 ) - (a: J + a> 2 ) (d+ y + z 2 )} + cz, + kz 2 = 0 


+ a 

cos 

0 [2qa> 1 - 

("2 - W 3 W l)] 



+ a 1 

sin 

B [2qa> 2 + (wj + w 2 a; 3 )] 



- O 

2 + o> 2 ) (d' + 

y + *1)} + c' 

+ k' z[ = 0 

(12) 

+ a' 

cos 

6> [2qcu 2 

+ + a> 2 a> 3 ) 

•] 


- a' 

sin 

0 [2qwj 

- (^ 2 - " 3 ^ 1 ) 

] 


- (<^x + 

co\) (d' + 

y + * 2 )} + c< 

zj + k' zj = 0 

(13) 


In trying to apply the Principle of Stability in the First Approximation, we 
linearize the system equations of motion in a neighborhood of the desired equi- 
librium point defined by 


co 3 r n , constant 


'l = W 2 = Z 1 = Z 2 = Z 1 = Z 2 = 0 


We substitute co, = n + v, and from the linearized version of (28), see that 
v 3 = 0. Without loss of generality, we put v 3 = 0. The other equations, when 
linearized, yield 

l i "l -{^2 - I 3 > n - J 3 <l} w 2 + ma ( z 2 + n2 

+ m'a' {cos £ (z 2 + (n +q) 2 z 2 ) + sin 0(z x ' + (n + q) 2 z|)} = 0 (15) 

l 2*>2 " {(*3 " V n + J3 fl}"! * ma (*! + n2 z j) 

-m'a' {cos 6(z’j + (n + q) 2 zj) - sin 0(z* 2 + (n + q) 2 z')} = 0 
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(17) 


m (y + *z'j - a(-v 2 - r.ojj)} + cZj + kZj = 0 

m{y + z’ 2 + a(Wj + nw 2 )| + cz 2 + kz 2 = 0 


(18) 


m' { y + z ' - a' cos 6 (<1> 2 - (n + 2q) Wj) 

+ a' + (n + 2q)w 2 )| + c'zj + k'zj = 0 (19) 

m' { y + z' + a' cos # + (n + 2q) w 2 ) 

+ a' sin £ (i- 2 - (n + 2q) a ’ 1 )} + c ' z' + k' z 2 = 0 (20) 


The stability of the origin of this set of equations has to be examined. 

A SIMILARITY TRANSFORMATION FOR LIAPUNOV REDUCTIBILITY 

Assume that the system of equations (15) -(20) is put in a standard vector- 
matrix first order differential equation form 

x = A(t) x (21) 

where x is an appropriately defined state vector in R ,0 and A(t) is the 10 x 10 
system matrix which has periodic coefficients. The Liapunov Reducibility 
Theorem 5 asserts that there exists a bounded nonsingular transformation P(t) 
with periodic coefficients with a bounded inverse P -1 (t) such that the matrix B 
defined by the relation 

B = P" 1 AP - P- 1 P (22) 

is a constant matrix. We may then define 

x=P(t)x (23) 
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and obtain 


x = Bx (24) 

The autonomous system defined by (24) is kinematically similar to the original 
system (21) and has similar stability behavior. 5 We may therefore confine our 
attention to the examination of the stability of (24). This is easy to do because 
the Routh-Hurwitz criterion provides a well defined algorithm. 

The handling of the tenth order system matrix defined in (21) can be clumsy 
and we would rather use the scalar equations (15)-(20) directly. To simplify the 
equations we introduce the approximation that y is identically zero. It can be 
seen from the definition of y in (1) that this is equivalent to assuming that the 
mass of any particular damper is negligible compared to the total mass of the 
system. This assumption implies that the location of the center of mass of the 
system is fixed in the body and is unaffected by the vibrations of the damper 
masses. We also define 

^n+Jjq-Ijn 


5 n + J3 q- I 2 n 
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ma 


ma . hi a T t . ma _ T ( 

— = L, ; — - ~ L 2» 

X 1 A 1 A 2 


m # a' 


m a 


_=L' 
A <■> 


c/m = c; c'/m' = c 1 ; k/m = k ; k'/m' = k' , 

n + q = r.j ; n + 2q = n 2 

rhe equations (15) -(20), may now be written 

+ \ 2 u 2 + Lj(z 2 + n 2 z 2 ) + L', {cos5(z' + n 2 z 2 ) 

+ s in 0 (z'j + n 2 z|)} = 0 


"a - A. jO)j - L 2 (z 1 + n 2 Zj) - L' 2 (cos^Czj + n 2 zj) 
- sin0(z 2 + n 2 z')} = 0 






L z d 



cos 6 sin 8 
-sin 8 cos 8 


L y 2j 


we have 


"l + + L i ( z 2 + n2 z 2 ) + L i (y 2 " 2c lyi + nn 2 y 2 ) = 0 

w 2 - \jCJj - L 2 (z\ - n 2 z j ) - L 2 (y'j + 2qy 2 + nn 2 y t ) = 0 
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( 25 ) 


(26) 


(27) 

(28) 
(29) 


Zj - a (a> 2 - nwj) + cz 2 + kZj = 0 

• A # A 

z 2 + a(co 1 + no> 2 ) + CZ 2 + kz 2 = 0 
z‘j - a' cos 6 (o> 2 - njCUj) + a' sin^(wj + n 2 o; 2 ) + c'zj + k'zj = 0 (30) 

*a “ a ’ cos ^ ("i + n 2 w 2 ) + a' sin d(u 2 - n^) + c ' z' + k’ z 2 = 0 (31) 

Introducing the nonholonomic transformation 

! r 
y t 


(32) 


(33) 

(34) 




iiiriiifiiifBiTiir- 


Zj - a (<i 2 - nc^j) + cij + k Zj =0 (35) 

z' 2 + a (ojj + nwj) + c z 2 + (< z 2 =0 (36) 

y‘l + 2q y 2 - q 2 y, - a* (w 2 - n 2 coj) + £' (y t + qy 2 ) + k' y l =0 (37) 

y 2 " 2qyi - q 2 y 2 + a ' ("i + n 2 ^ 2 ) + e' (y 2 - qy t ) + k' y 2 = 0 (38) 


(37) is obtained by multiplying (30) by cos # and (31) by -sin# and adding; (38) is 
similarly obtained by multiplying (30) by sin# and (31) by cos# and adding. 
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The transformation may clearly be put in the more elegant form suggested 
in (21). However, the above presentation is expected to bring out the basic idea 
more distinctly. The matrix in (32) simply defines a rotation operation and is 
extremely well known and widely used. In the present context, this transforma- 
tion effectively replaces the two rotor dampers by two equivalent ones which are, 
as it were, fixed in the body frame. The mutual coupling terms therefore have 
constant coefficients. 

AN EXAMPLE: THE SAS-A SPACECRAFT 

SAS-A is a dual-spin spacecraft launched at the end of 1970. It has been 
subjected to exhaustive study, a recent work being that of Bainum. 7 It has a 
pendulus nutation damper on the main body and a slightly flexible momentum 
wheel. Bainum takes axis 2 (Y axis) as the common spin axis and derives the 
nonlinear equations of motion. When linearized in a neighborhood of the assumed 
equilibrium motion, the following equations are obtained: 

B i>. + mr. (r. + r 0 )<6 =0 (41) 


Ao>j + fia> 3 (C - B) - co 3 I R 


s - 2m r j 



+ I Ri ( [a y + (fi + s) 2 a y ] cos 8 + [a^ + (fl + s) 2 a t ] sin 8} = 0 


(42) 
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Co> 3 + flojj (B - A) + u> 1 I R s - mr, <f> + m r. Q 2 4‘ 

2 M M 


t 

t 

+ I Rj (- (a y + (Q + s) 2 a y ] sin 9 + [a x + (0 + s) 2 a f ] 

cos (9} 

= 0 (43) 

K 

| 




4. 

I 

mr l ( 1 - - ) w - mr, — a i mr (r Q + r ) v + mr r 


I n 2 ^> 

| 

J 

\ M / M \ 

M J 


J 

’t 

| 

9 w , 

+ mr. + kd> + K0 = 0 


(44) 

i 

M 



i 

i 




* 

i 

i 

T Rl *, + k R ®, + [Kr + Ir, ( n + s) 2 ] a t 




+ { [o) 3 + (f2 + 2s) ojj J cos 9 + [cl'j - (H + 2s) i» 3 ] sin (9} = 0 (45) 

I R, % + k R% + 1 *r + *r, < n + s ) 2 ] % 

+ I Rj (- (w 3 + (fi + 2s) Wj] sin 5 + [cLj - (H + 2s)a.' 3 ] cos0} = 0 (4®) 

where Bainum's 7 notation has been used with the addition of 

o> 2 = Q + v 2 , fl constant 



Eq. (41) may be eliminated by solving for i> 2 and substituting in (44). 
Bainum introduces the following notations and approximations 

A = C = A 
B = B 

r = -tM/M 


k = {n(B - A) + I B s }/A 
m/M « 1 


We introduce a transformation 



cos 9 - s in 9 


a 

z 

s in 9 cos 9 


a 

L y_l 


(43) 


(49) 


The reduced set of kinematically similar autonomous equations are 


| m(r 0 + rj) 2 ]. .j k ■ f K r o fi2 1 

V — s — | 


</> 


r . r n A 

o>, + — l2a). = 0 

- o r 1 

r l r l 


(50) 


2mr,r . \A A . .. 

— fl<£ + - — cj 3 - - — (ij + 2 s /3j - {/? 2 + H (fi + 2 s) /3 2 ) =0 (51) 
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The study of this set of equations is a standard exercise in linear systems 
theory. In particular, the stability of t; origin for the nominal parameter 
values may be easily studied by Routh's method. In the above, fi and s repre- 
sent respectively the main spacecraft spin rate and the rotor relative spin rate. 
We replace these by the symbols n and q respectively to be consistent with the 
notation used before. These are the only two parameters in the above equations 
that are adjustable: the nominal values being 1/12 RPM and 2000 RPM respec- 
tively. Using the values given by Bainum 7 for the other parameters, and de- 
noting the Laplace Transform complex frequency variable by s, we have the 
characteristic matrix Q(s) of the system as 
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The characteristic polynomial |Q(s) | may be evaluated for various assumed 
values of n and q using a digital computer. Application of Routh’s criterion 
shows that the system is stable for a reasonably wide range of spacecraft and 
rotor spin rates around the nominal values. This result is hardly surprising: 
flight data has shown the system to be stable and a number of preflight simu- 
lation and experimental studies and approximate analyses had drawn the same 
conclusion. This, however, is the first analytic proof of stability under steady 
state nominal and off-nominal operating conditions and depends only on the 
validity of the mathematical model. 





CONCLUSIONS 


It has been shown that as required by the Floquet Theory and the Liapunov 
Reducibility Theorem , 5 a similarity transformation does exist which reduces 
the linear time-varying system of variational equations of a dual-spin system 
to a linear time-invariant system and thus makes a complete and rigorous anal- 
ysis possible. The analysis has been applied to a representative mathematical 
model of a dissipative dual-spin system and also to a recent model 7 of the SAS-A 
spacecraft. The method of analysis makes use of the fact that most dual-spin 
spacecrafts have symmetrical roto - ' dissipations in the two transverse axes. 

It is evident that more detailed analytical studies can be made by use of the 
proposed approach. Conclusive proofs for steady state spin stability may be 
obtained. Moreover, since the transformation is relatively easy to apply, the 
steady state spin stability for various assumed spin rates over a wide range of 
values may be studied. If the property of quasi-stationarity is assumed (frozen- 
time approximation), the above results may be used to analyze stability during 
slow h pin -up, a problem nut considered in any detail in present designs. 
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